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The authors attempt to establish the relative biological effectiveness (RBE) calculation for designing therapeut- 
ic proton beams on the basis of microdosimetry. The tissue-equivalent proportional counter (TEPC) was used 
to measure microdosimetric lineal energy spectra for proton beams at various depths in a water phantom. An 
RBE-weighted absorbed dose is defined as an absorbed dose multiplied by an RBE for cell death of human 
salivary gland (HSG) tumor cells in this study. The RBE values were calculated by a modified microdosi- 
metric kinetic model using the biological parameters for HSG tumor cells. The calculated RBE distributions 
showed a gradual increase to about 1cm short of a beam range and a steep increase around the beam range 
for both the mono-energetic and spread-out Bragg peak (SOBP) proton beams. The calculated RBE values 
were partially compared with a biological experiment in which the HSG tumor cells were irradiated by the 
SOBP beam except around the distal end. The RBE-weighted absorbed dose distribution for the SOBP beam 
was derived from the measured spectra for the mono-energetic beam by a mixing calculation, and it was 
confirmed that it agreed well with that directly derived from the microdosimetric spectra measured in the 
SOBP beam. The absorbed dose distributions to planarize the RBE-weighted absorbed dose were calculated 
in consideration of the RBE dependence on the prescribed absorbed dose and cellular radio-sensitivity. 
The results show that the microdosimetric measurement for the mono-energetic proton beam is also useful for 
designing RBE-weighted absorbed dose distributions for range-modulated proton beams. 
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INTRODUCTION 

Charged particle beams can deliver a relatively high dose 
to a target volume deeply seated in a patient's body better 
than conventional photon beams can, by means of their 
dose distribution along a Bragg curve [1]. Cancer therapy 
with high-energy proton beams was started at the Lawrence 
Berkeley National Laboratory in the United States in 1954. 
To date, proton beam therapy has been performed as a part 
of modern radiotherapy at about 30 facilities in many coun- 
tries, and the number of new proton facilities is rapidly in- 
creasing around the world. 



In particle radiation therapy, it is essential to calculate 
not only the absorbed dose but also the biological effect, 
which is often expressed by relative biological effectiveness 
(RBE). Most proton treatment facilities have adopted the 
constant RBE of 1.1 in clinical use [2, 3]. That means 
protons are assumed to be 10% more effective than refer- 
ence photons under any beam conditions, because the ac- 
curate measurement and calculation of the RBE have been 
difficult due to the complex dependence of biological 
effects on the dose level, radiation quality, biological end- 
point and cellular conditions. However, many biological 
experiments indicate that the RBE of low-energy protons 
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with high Hnear energy transfer (LET) is higher than that of 
high-energy protons with low LET [4, 5]. Paganetti et al. 
calculated the RBE for proton beams using the biological 
weighting function based on the microdosimetric lineal 
energy spectra [6]. Wilkens and Oelfke proposed a simple 
phenomenological RBE model for the fast calculation of 
proton RBE from dose-averaged LET in treatment planning 
[7]. Kundrat presented the probabilistic two-stage model 
with model parameters dependent on the particle type and 
LET value [8]. 

The microdosimetric approach is useful for assessing 
biological effects on the basis of the energy deposition pat- 
terns on microscopic targets [9], and it was applicable to es- 
timate the RBE for particle treatment beams [10-13]. The 
microdosimetric kinetic model (MKM) was developed by 
Hawkins to calculate cell survival curves for various types 
of radiation beam [14]. The MKM was then modified to 
derive the RBE for particle- and energy-modulated beams 
from the microdosimetric spectrum measured using the 
tissue-equivalent proportional counter (TEPC) [15]. The 
MKM was often used for calculating the RBE distributions 
in combination with the Monte Carlo calculation for 
photon beams [16], passive carbon ion beams [17, 18] and 
scanning carbon ion beams [19]. 

The purpose of this study was to derive an RBE- 
weighted absorbed dose distribution from microdosimetric 
spectra measured by the TEPC by means of the modified 
MKM. The radiation-induced cell death of human salivary 
gland (HSG) tumor cells was adopted as an endpoint to 
calculate the RBE, in accordance with the RBE model for 
carbon ion therapy at the National Institute of Radiological 
Sciences (NIRS) in Japan [20]. The RBE- weighted absorbed 
dose distribution of range-modulated proton beams was 
derived with the mixing calculation from the measured quan- 
tities for a mono-energetic proton beam before the modula- 
tion. The mixing calculation method was examined by 
comparison with the direct TEPC measurement of the range- 
modulated beam. Some depth-dose distributions with con- 
formal RBE-weighted absorbed doses in the spread-out Bragg 
peak (SOBP) regions were designed to demonstrate the de- 
pendence of RBE on the SOBP width, the prescribed dose 
and the tissue response (linear-quadratic a/p ratio for X-rays), 
in assuming that range modulation could be performed by the 
active beam scanning method or the layer-stacking method. 

MATERIALS AND METHODS 

Microdosimetric lineal energy measurement 

The microdosimetric spectra of proton beams were mea- 
sured using a proportional counter (LET- 1/2; Far West 
Technology, Inc., Goleta, CA, USA). The sensitive volume 
of the counter was nominally a sphere of 12.7 mm 
diameter, around which a spherical wall of A- 150 
tissue-equivalent plastic (thickness 1.27 mm) was situated. 



The wall was housed in an aluminum shell (thickness 
0.178 mm) to maintain a proportional gas under low pres- 
sure. The propane-based tissue-equivalent gas (p-TEG; 
54.6% C3H8, 40.16% CO2 and 5.26% N2, by volume) was 
enclosed with a pressure of 4.4 kPa in the counter to simu- 
late the energy imparted to a spherical tissue with a unit 
density of 1.0 jim diameter. The electric impulse from the 
counter entered a preamplifier (model 142B; ORTEC, Oak 
Ridge, TN, USA), and then the pre-amplified output signal 
was divided to enter three linear amplifiers (model 572/ 
671; ORTEC) with three gains of 2, 30 and 450 times to 
realize a wide measurement range and fine energy reso- 
lution. These signals were stored in multi-channel analyzers 
(MCA8000A; AMPTEK, Bedford, MA, USA). A calibra- 
tion pulser was used to examine the linearity between the 
pulse height of the pre-amplifier output and the channel 
and to relate those channels of the three spectra. The bin 
sizes of energy imparted were 3.0, 45 and 671 eV for 
high-, middle- and low-gain spectra, respectively. The 
lowest levels of energy imparted were 15 eV, 4.4 keV and 
67 keV for high-, middle- and low-gain spectra, respective- 
ly. Dead times were suppressed to less than 1%. 

Microdosimetry for proton beams 

The TEPC measurement was performed for a 155-MeV 
mono-energetic proton beam and a therapeutic proton beam 
with a nominal SOBP width of 6 cm. Figure 1 shows a 
schematic layout of the measurement in the proton treat- 
ment room with the rotating gantry in the Proton Medical 
Research Center of the University of Tsukuba. The gantry 
angle was set at 270° to use the horizontal beam in the 
TEPC measurement with a water phantom. The proton 
beam was laterally broadened by the double-scattering 
method before entering the treatment room. The field size 
was about 20 cm in diameter. A ridge filter composed of 
an array of metallic ridge bars was used to make the SOBP 
beam. Each ridge bar was designed to be a specific step- 
wise shape so that the absorbed dose was flattened in the 
SOBP region. 

Microdosimetric spectra were measured by the TEPC at 
various depths in the water phantom. The count rate of the 
TEPC was reduced to about lO"^ counts/s by reducing 
the beam intensity. The radiation quality of the low-intensity 
beam was assumed to be the same as that of the usual beam 
intensity. A main dose monitor of a parallel-plate ionization 
chamber was positioned in the beam line, but it was not 
available in the TEPC measurement due to the very low in- 
tensity of the beam. Thus, a small scintillation detector was 
placed on the edge of the water phantom to count an area of 
protons near the edge of the low-intensity beam. 

RBE calculation by the modified MKM 

The RBE-weighted absorbed dose was defined as a product 
of the absorbed dose and RBE for cell death of in vitro 
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Fig. 1. Schematic layout of the irradiation system for the therapeutic proton beam and TEPC measurement. 



HSG tumor cells. The HSG tumor cells have been used 
to determine the RBE in carbon ion beams as the standard 
reference cell line [20, 21] and to research the RBE in 
proton and carbon ion beams [22-24]. The RBE value was 
calculated from dose-survival curves of radiation of interest 
and reference photon radiation. The surviving fraction, S, 
was calculated from the lineal energy spectrum by the 
modified MKM with the model parameters of the HSG 
tumor cells (ao = 0A3Gy~\ (3 = 0.05 Gy"^ rj = 0.42|im 
and yo= 150 keV/|im) as follows [15, 25]: 



HSG tumor cells is defined as follows: 



S = exp 



/3 



D- j5D^ 



(1) 



RBE-D = ^^^^^^^^ (4) 
2/3 ' V y 



where the unit of Gy (RBE) was used for the D^be value. 

The model parameters given above were assumed to be 
radiation-independent values for HSG tumor cells, and ori- 
ginally obtained to calculate RBE values for HSG tumor 
cells from microdosimetric spectra measured by the TEPC. 
Meanwhile, the different model parameters were used in 
other papers [19, 26] since those were adjusted to calculate 
survival curves from averaged energies imparted to the cy- 
lindrical volume using amorphous track structure models. 



ygj(l-exp(-//yg))f(y)dy 



(2) 



where D is the absorbed dose, p is the density of tissue 
assumed to be p =1 g/cm^, f{y) is the probability density 
of lineal energy, y, represents the saturation-corrected 
dose-mean lineal energy and P is the constant value of 
0.05 Gy-^. 

The RBE value of HSG tumor cells with S was calcu- 
lated as follows: 



RBE 



Ds,R 
Ds 



^/al-4|3ln{S)-aR 
2f5Ds 



(3) 



where Ds and Ds^r are the doses necessary for S by radi- 
ation of interest and reference radiation, respectively. For 
this RBE calculation, the 200-kV X-rays were used as the 
reference radiation, which had the a value of 0.19 Gy~^ 
and the (3 value of 0.05 Gy~^ in the linear quadratic model 
[15]. Then, the RBE- weighted absorbed dose, Drbe^ for 



Mixing calculation for range-modulated beams 

The therapeutic SOBP beam is made by a superposition of 
many range-modulated Bragg curves. The mixed absorbed 
dose, Dfn, and mixed value, j*^, of the range-modulated 
beam were calculated from those of the mono-energetic 
beams by the following mixing calculation [25]: 



Dfn{x) 



^y^{x^ Si)di{x^ Si)ri 



y^m{x) — 



^di{x^Si)ri 



(5) 



(6) 



where di(x) and are the absorbed dose and value, 
respectively, at a depth of x in the mono-energetic beam, 
and is the weighting factor required for the superimpos- 
ition with the shifted depth of Sj. The survival fraction and 
the RBE-weighted absorbed dose in the range-modulated 
beam were calculated from the absorbed dose and distri- 
butions in the mono-energetic beam by Eqns (l)-(6). 



488 



Y. Kase et al 



Cell irradiation 

The HSG tumor cells (JCRB1070: HSGc-C5) were used in 
the measurement of survival curves. The cells were seeded 
to about 5 X 10^ cells per 25 -cm^ culture bottle (Nunc, 
Roskilde, Denmark) at the time of irradiation. Eagle's 
minimum essential medium (E-MEM) supplemented with 
10% fetal bovine serum (FBS; JRH Biosciences, Lenexa, 
KS, USA) and antibiotics (penicillin and streptomycin) was 
used at cell culture. Seven bottles were prepared to obtain 
one survival curve. One bottle was not irradiated, to ascer- 
tain the plating efficiency of the series of colony formation 
experiments, and the other bottles were irradiated once by 
the therapeutic proton beam with dilferent absorbed doses. 
Survival curves of the HSG tumor cells were measured at 
four depths of entrance, 5 mm in the back of the proximal 
SOBP peak, at the center of the SOBP and 5 mm short of 
the distal SOBP peak in the therapeutic proton beam. 

After being exposed to radiation, the cells in the irra- 
diated bottles were harvested by trypsinization and counted 
using a cell counter (Coulter, Z-1; Beckman Coulter Inc., 
Brea, CA, USA). Then the cells were diluted with E-MEM 
and plated onto three dishes of 6-cm diameter (Falcon, 
3002; BD, Franklin Lakes, NJ, USA) to yield about 100 
colonies per dish. After 2 weeks of incubation, the colonies 
were fixed and stained. Colonies containing more than 50 
cells were counted as surviving colonies. The average 
plating efficiency of control experiments was about 81%. 

RESULTS 
Microdosimetric spectra 

The microdosimetric lineal energy spectrum, yd{y), in the 
semi-log representation is the standard representation of a 
microdosimetric spectrum. Figure 2 shows typical yd{y) 
spectra measured by the TEPC at typical depths in the 



water phantom for 155-MeV proton beams of the mono- 
energetic beam and the typical treatment beam with the 
6-cm nominal SOBP width. The lineal energies of the 
proton beam ranged from 0. 1 to 20 keV/|im around the en- 
trance, and from 1 to 100 keV/jim around the Bragg peak. 
The ydiy) spectra shifted to higher y events with the depth 
up to the maximum beam range because of the decreased 
proton energy and the increased average LET of the proton. 

Figure 3 shows the frequency-mean lineal energy, yp, 
dose-mean lineal energy, y^), and saturation-corrected dose- 
mean lineal energy, y*, with the saturation parameter 
y 0=150 keV/|im as a function of the depth in water, mea- 
sured by the TEPC for the mono-energetic 155-MeV 
proton beam. The y^ values measured by the TEPC often 
have large uncertainty because they are significantly 
alfected by even a few events with large lineal energy. The 
high-energy proton sometimes undergoes a nuclear reaction 
in material. The nuclear reaction may produce recoil neu- 
trons and radionuclides, which can cause enormously large 
lineal energy. On the other hand, the y* values were rela- 
tively stable results, because of oversaturation with the 
elfect of lineal energy larger than 150 keV/|im along with 
the radiobiological overkill effect. The events beyond the 
proton maximum range are also attributed to the neutrons 
and gamma-rays emitted from the radionuclides [13]. 

RBE-weighted absorbed dose distributions 

The depth-dose and depth-y* distributions were calculated 
from the microdosimetric spectra measured by the TEPC. 
Figure 4 shows the depth-dose distributions and depth- 
RBE distribution for the mono-energetic 155-MeV proton 
beam. The values were around 4 keV/|im in the dose 
plateau region from the entrance to the depth of 1 cm short 
of the dose peak, and they increased sharply to about 18 
keV/|im beyond 6 mm after the dose peak. The RBE values 
obtained by the TEPC measurement were almost constant 
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Fig. 2. Microdosimetric spectra of lineal energy, yd{y), measured using the TEPC at typical depths in the water phantom for the 
155-MeV proton beams: (a) mono-energetic beam; (b) typical treatment beam with a 6-cm nominal SOBP width. 
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from the entrance to the depth of 1 cm short of the dose 
peak, and they increased sharply with the depth in water 
around the dose peak. The calculated RBE values were 
0.95 ±0.02 at the entrance (22 mm), 1.05 ±0.01 at the 
proximal 90% dose depth (128 nmi), 1.16 ±0.01 at the 
dose peak (131 mm), 1.22 ±0.01 at the distal 90% dose 
depth (132 mm), 1.30 ±0.02 at the 50% dose depth (133 
mm) and 1.74 ± 0.07 at the 10% dose depth (136 mm). The 
uncertainties in the RBE value were calculated from the un- 
certainties in the statistical uncertainties of the absorbed 
dose and the value. The calculated RBE values were 
below 1.0 in the plateau region up to about 100 mm since 
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Fig. 3. Frequency-mean lineal energy, yp, dose-mean lineal 
energy, y^, and saturation-corrected dose-mean lineal energy with 
the saturation parameter jo= 150 keV/|im, j*, as a function of the 
depth in water, measured by the TEPC for the mono-energetic 
155-MeV proton beam. 



the values measured in the region were less than about 
4.3 keV/|im corresponding to reference X-rays [15]. 

Figure 5 shows survival curves of the HSG tumor cells 
for the therapeutic proton beams with 6-cm nominal SOBP 
width. The survival curve for the 200-kV X-rays [16] was 
also plotted to derive experimental RBE values for the 
proton beams. If these survival curves were fitted by the 
linear quadratic model, then the fitted a and P values were, 
respectively, 0.157 Gy~^ and 0.048 Gy~^ at the entrance, 
0.164 Gy"^ and 0.048 Gy"^ at 5 mm in the back of the 
proximal SOBP peak, 0.139 Gy"^ and 0.056 Gy"^ at the 
center of the SOBP, 0.141 Gy"^ and 0.062 Gy"^ at 5 mm 
short of the distal SOBP peak in the therapeutic proton 
beam and 0.196 Gy"^ and 0.044 Gy"^ for 200-kV X-rays. 
However, fitting two parameters with these limited survival 
data would cause large errors because of the anticorrelation 
between a and P values. Thus, only a values were fitted 
with a constant P value in analogy with the MKM calcula- 
tion to calculate RBE values. The a values fitted with the 
constant p value of 0.05 Gy"^ were 0.147 ±0.005 Gy"^ at 
the entrance, 0.155 ± 0.004 Gy~^ at 5 mm at the back of the 
proximal SOBP peak, 0.167 ±0.009 Gy"^ at the center of 
the SOBP and 0.194 ±0.007 Gy"^ at 5 mm short of the 
distal SOBP peak in the therapeutic proton beam and 
0.164 ±0.008 Gy"^ for 200-kV X-rays. 

Figure 6 shows depth-y* distributions for the 155-MeV 
proton beam with the 6-cm nominal SOBP width. Figure 7 
shows depth-dose distributions and the depth-RBE distri- 
bution for the same beam. The distributions calculated by 
the mixing calculation with Eqn (6) from the depth-y* dis- 
tribution for the mono-energetic beam were compared with 
those directly derived by Eqn (2) from the microdosimetric 
spectra measured in the SOBP beam. The directly derived 
RBE values agreed well with the calculated RBE distribu- 
tion within the uncertainties, except for the depth beyond 
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Fig. 4. Depth-dose distributions and depth-RBE distribution for the mono-energetic 155-MeV proton beam. Plots represent the TEPC 
measurement. The dashed line shows the measurement by a plane-parallel ionization chamber. The right figure shows an enlarged view 
around the beam range of the left figure. 
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Fig. 5. Survival curves of the HSG tumor cells for 155-MeV 
proton beams with the 6-cm nominal SOBP width and 200-kV 
X-rays. The cells were irradiated at four depths of entrance, 5 mm 
in the back of the proximal SOBP peak, at the center of the 
SOBP, and 5 mm short of the distal SOBP peak in the therapeutic 
proton beam. Polyethylene blocks were used for changing the 
water-equivalent depth. The lines represent the curves fitted by the 
liner quadratic model with the constant p value of 0.05 Gy~^. 
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Fig. 7. Depth-dose distributions and depth-RBE distributions 
for the 155-MeV proton beam with the 6-cm nominal SOBP 
width. The closed triangles, diamonds and squares represent the 
absorbed dose, the RBE-weighted absorbed dose and the RBE, 
respectively, directly derived by Eqn (2) from the microdosimetric 
spectra measured in the SOBP beam. The open squares represent 
the RBE obtained by the cell irradiation of the HSG tumor cells. 
The solid lines show the distributions calculated by the mixing 
calculation with Eqn (6) from the depth-dose and depth-y* 
distributions measured in the mono-energetic beam. 
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Fig. 6. Depth-y* distribution for the 155-MeV proton beam 
with the 6-cm nominal SOBP width. The closed circles represent 
the y* values directly measured in the SOBP beam. The solid 
lines show the distribution calculated by the mixing calculation 
with Eqn (6) from the depth-dose and depth-y* distributions 
measured in the mono-energetic beam. 



the maximum proton range. The calculated RBE values 
were 0.97 ± 0.02 at the entrance (22 mm), 1.00 ± 0.03 at the 
SOBP proximal region (73 mm), 1.03 ±0.03 at the SOBP 



center (98 mm), 1.09 ±0.02 at the SOBP rear region (123 
mm) and 1.23 ±0.09 at the beam range of the 90% dose 
(128 mm). The experimental RBE values were 0.95 ±0.05 
at the entrance, 0.98 ± 0.05 at the SOBP proximal region, 
1.01 ±0.05 at the SOBP center and 1.08 ±0.06 at the 
SOBP rear region 

Typically, the reproducibility of the experimental sur- 
vival fraction was poorer than that of the physical dose 
measurements. Some possible reasons for changing the 
survival fraction are unexpected changes of the radio- 
sensitivity of the HSG cells, the compositions of nutrient 
medium, culture conditions, and so on. It was assumed that 
the RBE value was a relatively reproducible result even 
though the reference survival curves changed, within 
about ± 5% from our previous experience. 

Designed depth-dose distributions 

Figure 8 shows depth-dose distributions for 155-MeV 
proton beams with the SOBP widths of 2, 4, 6, 8 and 10 
cm, which were designed to be constant RBE-weighted 
absorbed doses of 2 Gy (RBE) in each SOBP region. The 
RBE increased with the depth and was almost dependent 
on the residual range in the SOBP region at any SOBP 
width. Therefore, to have a flat RBE-weighted absorbed 
dose distribution, the absorbed dose must decrease with the 
depth in the SOBP region. It is well known that the RBE is 
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Fig. 8. Depth-dose distributions for 155-MeV proton beams 
with the SOBP widths of 2, 4, 6, 8 and 10 cm, designed to be a 
constant RBE-weighted absorbed dose of 2 Gy (RBE) at the 
center of the SOBR The solid and dashed lines represent the 
absorbed dose and RBE-weighted absorbed dose, respectively. 
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Fig. 9. Depth-dose distributions for 155-MeV proton beams 
with a 6-cm SOBP width for the cases of RBE-weighted absorbed 
doses of 0.5, 1, 2, 5 or 10 Gy (RBE), designed to be a constant 
RBE-weighted absorbed dose in the SOBP region. The solid and 
dashed lines represent the absorbed dose and RBE-weighted 
absorbed dose, respectively. 



dependent on the absorbed dose and radio-sensitivity for 
the irradiation target. Figure 9 shows calculated depth-dose 
distributions for the 155-MeV proton beams with a 6-cm 
SOBP width for various prescribed doses, which were 
designed to be constant RBE-weighted absorbed doses of 
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Fig. 10. Depth-dose distributions for 155-MeV proton beams 
with a 6-cm SOBP width for the cases of alp ratios of 1, 2, 5, 10 
and 20 Gy for X-rays with the constant p value of 0.05 Gy~^, 
designed to be a constant RBE-weighted absorbed dose of 2 Gy 
(RBE) in the SOBP region. The solid and dashed lines represent 
the absorbed dose and RBE-weighted absorbed dose, respectively. 



0.5, 1, 2, 5 or 10 Gy (RBE) in the SOBP region. Figure 10 
shows depth-dose distributions for the 155-MeV proton 
beams with a 6-cm SOBP width for various a/p ratios of 1, 
2, 5, 10 and 20 Gy for 200-kV X-rays with the constant (3 
value of 0.05 Gy~^, which were calculated in the same 
manner as in a previous study [25]. They were designed to 
be constant RBE-weighted absorbed doses of 2 Gy (RBE) 
in the SOBP region. 

The larger the absorbed dose or the larger the a/p ratio 
for X-rays was, the lower the RBE derived by the modified 
MKM calculation, notably in the SOBP region. Therefore, 
the absorbed dose distribution needs to be optimized 
depending on the prescribed absorbed dose and cellular - 
to result in a flat RBE-weighted absorbed dose distribution 
in the SOBP region. 

DISCUSSION 

The RBE distributions calculated from the TEPC measure- 
ment by the modified MKM showed a gradual increase to 
about 1 cm short of a beam range and a steep increase 
around the beam range for both the mono-energetic and 
SOBP proton beams. These results were roughly consistent 
with the calculated RBE distributions reported in several 
previous papers [6, 7, 27]. Matsuura et ah [28] reported that 
the experimental RBE value of HSG tumor cells at the 
Bragg peak for the mono-energetic 235-MeV proton beam 
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was 1.18 ± 0.07 times higher than that at the plateau region 
and nearly independent of the dose rate. The extension of 
the beam range by about 1 mm was also observed, match- 
ing what Carabe et al reported [29]. However, a systemati- 
cal uncertainty in the water-equivalent depth of about 
±1.5 mm was estimated in consideration of the spherical 
wall and placement accuracy of the TEPC in the water 
phantom. A miniature TEPC would be effective for measur- 
ing microdosimetric spectra around the beam range with 
higher resolution [27]. 

We validated the calculation method to derive the 
RBE-weighted absorbed dose distribution of the range- 
modulated proton beam from the microdosimetric spectra 
of the mono-energetic proton beam with the initial energy 
of 155 MeV. The result indicates that the microdosimetric 
measurement is available not only to evaluate the 
RBE-weighted dose distribution for designed SOB? beams, 
but also to design new SOBP beams from elemental mono- 
energetic proton beams. 

The mixing calculation in this study assumed that the 
value and absorbed dose distribution of the same residual 
range were approximately unchanged by the range shift. 
Meanwhile, the RBE distribution is known to depend on 
the initial beam energy. The larger the initial beam energy 
was, the smaller the maximum RBE around the beam range 
was [6]. Therefore, if different initial beam energies are 
used for designing the RBE-weighted dose distribution, the 
value and absorbed dose distribution should be mea- 
sured or calculated with respect to each initial beam energy. 

The change tendency of the RBE distribution depend- 
ence on the physical radiation quality is thought to be ap- 
proximately common to that of other human cell lines. The 
variable RBE calculated by the modified MKM has been 
applied to scanning carbon ion treatment planning [19]. In 
the future, the variable RBE calculation will be important 
to achieve more high-precision proton therapy. An active 
beam scanning method and layer- stacking method would 
make it possible to irradiate treatment beams with consist- 
ently flat RBE-weighted absorbed dose distribution in bio- 
logically modulated particle therapy. 

Incidentally, the absolute RBE value obtained for the 
in vitro HSG tumor cells should not be simply applied in 
clinical situations, because it would be greatly influenced 
by cellular characteristics, circumstance, culture solution, 
and so on. Therefore, it is difficult in this study to discuss 
an absolute dose prescription or a clinical RBE for proton 
beam therapy. Moreover, the absolute RBE value is also 
dependent on the choice of reference radiation type. The 
survival curve of a photon beam depends on the photon 
energy spectrum. For example, the absorbed doses of 10% 
survival level for the HSG tumor cells were 5.3 Gy for 
200-kV X-rays, 5.9 Gy for ^^Co gamma-rays, and 6.1 Gy 
for 6-MV X-rays [30]. Therefore, the RBE values calcu- 
lated at a 10% survival level usins the '^°Co gamma-rays 



and 6-MV X-rays as reference radiation are about 1 1 % and 
15% larger, respectively, than that using the 200-kV 
X-rays. Therefore, this study result should be used for im- 
proving only relative dose distribution in a proton beam 
treatment planning system. 

In conclusion, we calculated the RBE of HSG tumor 
cells by the modified MKM using microdosimetric spectra 
measured by the TEPC for 155-MeV proton beams. The 
calculated RBE values were partially compared with the 
biological experiment with HSG tumor cells irradiated by 
the proton SOBP beam except around the distal peak, 
though there was the causeless shift of the survival curve 
for reference X-rays. The mixing calculation was effective 
to derive the RBE-weighted absorbed dose distribution of 
the range-modulated proton beam from the microdosimetric 
spectra of the mono-energetic proton beam with the same 
initial proton energy. The modified-MKM calculation 
showed a sharp increase of the RBE of about ±1 cm 
around the proton beam range, though the RBE values 
around the dose peak were not confirmed by cell irradi- 
ation. The absorbed dose distributions to planarize the 
RBE-weighted absorbed dose were calculated in consider- 
ation of the RBE dependence on the prescribed absorbed 
dose and cellular radio-sensitivity. We expect that this 
method would make it possible to calculate the 
RBE-weighted dose distributions easily with the microdosi- 
metric quantities for various irradiation conditions in a 
future treatment planning system. 
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